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QCD, Gauge Fixing, and the Gribov Problem 
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The standard techniques of gauge-fixing, such as covariant gauge fixing, are entirely adequate for the purposes 
of studies of perturbative QCD. However, they fail in the nonperturbative regime due to the presence of Gribov 
copies. These copies arise because standard local gauge fixing methods do not completely fix the gauge. Known 
Gribov-copy-free gauges, such as Laplacian gauge, are manifestly non-local. These issues are examined and the 
implications of non-local gauge-fixing for ghost fields, BRST invariance, and the proof of renormalizability of 
QCD are considered. 
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1. GAUGE FIXING 

The naive Lagrangian density of QCD is 

C QC B = -\F^F^ + J2<lf(^-mf)qf, (1) 
/ 

where the index / corresponds to the quark 
flavours. The naive Lagrangian is neither gauge- 
fixed nor renormalized, however it is invariant un- 
der local SU(3) C gauge transformations g(x). For 
arbitrary, small Lu a (x) we have 

g(x) = exp l-ig, (j) w a {x) J G SU{3), (2) 

where the A a /2 = t a are the generators of the 
gauge group SU (3) and the index a runs over the 
eight generator labels a = 1, 2, 8. 

Consider some gauge-invariant Green's func- 
tion (for the time being we shall concern ourselves 
only with gluons) 



(0| T(6[A]) |fi> = 



JVA Q[A ] e 
—fVA 



iS[A] 



,iS[A] 



(3) 



where 0[A] is some gauge- independent quantity 
depending on the gauge field, A IJi (x). We see that 
the gauge-independence of 0[A] and S[A] gives 
rise to an infinite quantity in both the numerator 
and denominator, which must be eliminated by 
gauge-fixing. The Minkowski-space Green's func- 
tions are defined as the Wick-rotated versions of 
the Euclidean ones. 
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Figure 1. Illustration of the gauge orbit contain- 
ing A^ and indicating the effect of acting on A^ 
with the gauge transformation g. The action S[A] 
is constant around the orbit. 



The gauge orbit for some configuration A^ is 
defined to be the set of all gauge-equivalent con- 
figurations. Each point A^ on the gauge orbit 
is obtained by acting upon A^ with the gauge 
transformation g. By definition the action, S[A], 
is gauge invariant and so all configurations on the 
gauge orbit have the same action, e.g., see the il- 
lustration in Fig. [|. 

The integral over the gauge fields can be writ- 
ten as the integral over a full set of gauge- 
inequivalent (i.e., gauge-fixed) configurations, 
jT>A g l , and an integral over the gauge group 
J T>g. In other words, J T>A s l ' is an integral over 
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the set of all possible gauge orbits and J T>g is an 
integral around the gauge orbits. Thus we can 
write 



VA 



VA 



g.f. 



Vg. 



(4) 



To make integrals such as those in the numera- 
tor and denominator of Eq. (^) finite and also to 
study gauge-dependent quantities in a meaning- 
ful way, we need to eliminate this integral around 
the gauge orbit, J T>g. 

2. GRIBOV COPIES AND THE 
FADDEEV-POPOV DETERMINANT 

Any gauge- fixing procedure defines a surface in 
gauge-field configuration space. Fig. || is a depic- 
tion of these surfaces represented as dashed lines 
intersecting the gauge orbits within this configu- 
ration space. Of course, in general, the gauge or- 
bits are hypersurfaces and so are the gauge-fixing 
surfaces. Any gauge- fixing surface must, by defi- 
nition, only intersect the gauge orbits at distinct 
isolated points in configuration space. For this 
reason, it is sufficient to use lines for the sim- 
ple illustration of the concepts here. An ideal 
(or complete) gauge-fixing condition, F[A] = 0, 
defines a surface that intersects each gauge or- 
bit once and only once and by convention con- 
tains the trivial configuration A^ — 0. A non- 
ideal gauge- fixing condition, F'[A] = 0, defines a 
surface or surfaces which intersect the gauge or- 
bit more than once. These multiple intersections 
of the non-ideal gauge fixing surface(s) with the 
gauge orbit are referred to as Gribov copies |Q- 
|(| . Lorentz gauge (d ll A tJ '(x) = 0) for example, 
has many Gribov copies per gauge orbit. By def- 
inition an ideal gauge fixing is free from Gribov 
copies. We refer to the ideal gauge-fixing surface 
F[A] — as the Fundamental Modular Region 
(FMR) for that gauge choice. Typically the gauge 
fixing condition depends on a space-time coordi- 
nate, (e.g., Lorentz gauge, axial gauge, etc.), and 
so we write the gauge fixing condition more gen- 
erally as F([A];x) = 0. 

Let us denote one arbitrary gauge configuration 
per gauge orbit as A^ and let this correspond to 
the "origin" of gauge configurations on that gauge 
orbit, i.e., to g = on that orbit. Then each 





F'[A]=0 



Figure 2. Ideal, F[A], and non-ideal, F'[A], 
gauge-fixing. 



gauge orbit can be labelled by A ^ and the set of 
all such A^ is equivalent to one particular, com- 
plete specification of the gauge. Under a gauge 
transformation, <?, we move from the origin of the 
gauge orbit to the configuration, A 9 , where by 

definition A° JU A% = gAtf - (i/g s )(d„g)gl 
Let us denote for each gauge orbit the gauge 
transformation, g = g[A°], as the transformation 
which takes us from the origin of that orbit, A^ , 
to the configuration, A&- 1 • , which lies on the ideal 
gauge-fixed surface specified by F([A];x) = 0. 
In other words, we have F([A];a;)|^a = for 
A3 = Af; L e FMR. 

The inverse Faddeev-Popov determinant is de- 
fined as the integral over the gauge group of the 
gauge-fixing condition, i.e., 



A^H* f -] = JvgS[F[A}} 

(5F{[A]-x) 



Vg S(g - g) 



det 



V 8g(y) 
Let us define the matrix Mj?[A] as 
SF a ([A};x) 



M F {[A];x,y) ab 



Sg b (y) 



(5) 



(6) 



Then the Faddeev-Popov determinant for an arbi- 
trary configuration A^ can be defined as Ap [A] = 
| det (The reason for the name is now 

clear). Note that we have consistency, since 
A^[^- f -] ee A^A*] = jVg 5(g-~g)Ap 1 [A}. 



3 



We have 1 = JVg A F [A] 5[F[A]] by definition 
and hence 



VA S 



.{. 



VA si J Vg A F [A] S[F[A}} 
J VA A F [A] 5[F[A}} 



(7) 



Since for an ideal gauge-fixing there is one 
and only one g per gauge orbit, such that 
F(L4];x)|g = 0, then |detM F L4]| is non-zero on 
the FMR. It follows that if there is at least one 
smooth path between any two configurations in 
the FMR and since the determinant cannot be 
zero on the FMR, then it cannot change sign on 
the FMR. The Gribov horizon is defined to by 
those configurations with detMpfA] = which 
lie closest to the FMR. By definition the determi- 
nant can change sign on or outside this horizon. 
Clearly, the FMR is contained within the Gribov 
horizon and for an ideal gauge fixing, since the 
sign of the determinant cannot change, we can 
replace jdetMp with det Mp, [i.e., the overall 
sign of the functional integral is normalized away 
inEq. (§)]. 

These results are generalizations of results from 
ordinary calculus, where 



det 



dfi 



dx 1 ---dx n 5 {n) {f{x)) (8) 



and if there is one and only one x which is a solu- 
tion of f(x) = then the matrix My = dfi/dxj 
is invertible (i.e., non-singular) on the hypersur- 
face f(x) = and hence det 



3. GENERALIZED 
TECHNIQUE 



FADDEEV-POPOV 



Let us now assume that we have a family of 
ideal gauge fixings F(L4];x) = f([A};x) — c(x) 
for any Lorentz scalar c(x) and for /([A]; x) be- 
ing some Lorentz scalar function, (e.g., d^A^x) 
or n^A^x) or similar or any nonlocal generaliza- 
tions of these). Therefore, using the fact that we 
remain in the FMR and can drop the modulus on 
the determinant, we have 



2M gi 



VA det M F [A\ S[f[A] 



c . 



(9) 



Since c(x) is an arbitrary function, we can define 
a new "gauge" as the Gaussian weighted average 
over c(x), i.e., 



PA g f - oc / Vc exp <J / d 4 xc{x^ 2 



x / VA det M F [A] 5[f[A] - c] 
oc J VA detMp[A] 



' / d 4 xf([A];x^ 



oc J VAV X Vx exp|-i J d 4 xd 4 y 
x x{x)M F ([A];x,y)x{y)} 
xexp^^Jd 4 xf([A];x)^ (10) 

where we have introduced the anti-commuting 
ghost fields x and x- Note that this kind of ideal 
gauge fixing does not choose just one gauge con- 
figuration on the gauge orbit, but rather is some 
Gaussian weighted average over gauge fields on 
the gauge orbit. We then obtain 



(n\T(d[...])\n) = 



J VqVqVAVxVx 0[.. 
j VqVqVAVxVx e iS ^ 



, (11) 



where 



S^[q,q,A, X ,x] = J d * 



4 ^ 



Z ^ f 



+ j d 4 xd 4 y X (x)M F ([A];x,y)x(y). (12) 

4. STANDARD GAUGE FIXING 

We can now recover standard gauge fixing 
schemes as special cases of this generalized form. 
First consider standard covariant gauge, which we 
obtain by taking /(L4];a;) = d fl A IJ, (x) and by ne- 
glecting the fact that this leads to Gribov copies. 
We need to evaluate M F [A] in the vicinity of the 
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gauge-fixing surface for this choice: 



M F ([A];x,y) ab = 



SF a ([A]:x) 

5g b {y) 



S[d^A a ^(x) - c(x)} 



d: 



.5A a ^{x) 



(13) 



5g»(y) " Sg»(y)] 

Under an infinitesimal gauge transformation g we 
have 

A *(x) sm ^ 9 (A%(x)=A;(x) 

+g s f abc Lo b {x)A c ll (x) - d^ a {x) + 0(lu 2 ) (14) 

and hence in the neighbourhood of the gauge fix- 
ing surface (i.e., for small fluctations along the 
gauge orbit around A&f-), we have 



M F ([A];x,y) 



ab 



9, 



JA^(x) 



(15) 



j=0 



= ( [-d x > J -S ab + g s f abc A c ^(x)] x 5^\x - y) 

We then recover the standard covariant gauge- 
fixed form of the QCD action 



Sc[q,q,A,x,x] 



d 4 x 



"4 \iu 



~{d^(x)f 



f 



q f (ip-m f )q f 



-{d^xa){d»5 ab - gf abc AZ) Xb ■ 



(16) 



However, this gauge fixing has not removed the 
Gribov copies and so the formal manipulations 
which lead to this action are not valid. This 
Lorentz covariant set of naive gauges corresponds 
to a Gaussian weighted average over generalized 
Lorentz gauges, where the gauge parameter £ 
is the width of the Gaussian distribution over 
the configurations on the gauge orbit. Setting 
£ = we see that the width vanishes and we ob- 
tain Landau gauge (equivalent to Lorentz gauge, 
d fJ 'A fl (x) = 0). Choosing £ = 1 is referred to as 
"Feynman gauge" and so on. 

We can similarly recover the standard QCD ac- 
tion for axial gauge, where n tl A ti (x) — 0. Pro- 
ceeding as for the generalized covariant gauge, we 
first identify = n^A^(x) and obtain the 

gauge-fixed action 



Si[q,q,A] 



d 4 x 



4 



~ (n^(x)f + q f (ip- mj)q f 



f 



(17) 



Taking the "Landau-like" zero-width limit £ — > 
we select n^A^{x) = exactly and recover the 
usual axial-gauge fixed QCD action. Axial gauge 
does not involve ghost fields, since in this case 



M F ([A^-];x,y) 



\ ab 



SA a ^(x) 



$co b (y)] 
-d x »8 ab ]5 {i \x-y) 



(18) 



which is independent of A^ since n^AZ (x) = 0. 
In other words, the gauge field does not appear 
in Mf[j4] on the gauge-fixed surface. Unfortu- 
nately axial gauge suffers from singularities which 
lead to significant difficulties when trying to de- 
fine perturbation theory beyond one loop. A re- 
lated feature is that axial gauge is not a complete 
gauge fixing prescription. While there are com- 
plete versions of axial gauge on the lattice, these 
always involve a nonlocal element, or reintroduce 
Gribov copies at the boundary so as not to de- 
stroy the Polyakov loop. 

5. DISCUSSION AND CONCLUSIONS 

There is no known Gribov-copy-free gauge fix- 
ing which is a local function of A fi (x). In other 
words, such a gauge fixing cannot be expressed 
as a function of A^(x) and a finite number of 
its derivatives, i.e., F(L4];x) ^ F(d ll ,A ll (x)) for 
all x. Hence, the gauge-fixed action, £$[•••], in 
Eq. ( |l2| ) becomes non-local and gives rise to a 
nonlocal quantum field theory. Since this action 
serves as the basis for the proof of the renormaliz- 
ability of QCD, the proof of asymptotic freedom, 
local BRS symmetry, and the Schwinger-Dyson 
equations (to name but a few) the nonlocality of 
the action leaves us without a reliable basis from 
which to prove these features of QCD in the non- 
perturbative context. 

It is well-established that QCD is asymptoti- 
cally free, i.e., it is weak-coupling at large mo- 
menta. In the weak coupling limit the func- 
tional integral is dominated by small action con- 
figurations. As a consequence, momentum-space 
Green's functions at large momenta will receive 
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their dominant contributions in the path integral 
from configurations near the trivial gauge orbit, 
i.e., the orbit containing = 0, since this orbit 
minimizes the action. If we use standard gauge 
fixing, which neglects the fact that Gribov copies 
are present, then at large momenta J VA will be 
dominated by configurations lying on the gauge- 
fixed surfaces in the neighbourhood of each of 
the Gribov copies on the trivial orbit. Since for 
small field fluctuations the Gribov copies cannot 
be aware of each other, we merely overcount the 
contribution by a factor equal to the number of 
copies on the trivial orbit. This overcounting is 
normalized away by the ratio in Eq. (||) and be- 
comes irrelevant. Thus it is possible to under- 
stand why Gribov copies can be neglected at large 
momenta and why it is sufficient to use standard 
gauge fixing schemes as the basis for calculations 
in perturbative QCD. Since renormalizability is 
an ultraviolet issue, there is no question about 
the renormalizability of QCD. 

Lattice QCD has provided numerical confirma- 
tion of asymptotic freedom, so let us now turn 
our attention to the matter of Gribov copies in 
lattice QCD. Since the observable 0[A] and the 
action are both gauge-invariant it doesn't mat- 
ter whether we sample from the FMR of an ideal 
gauge-fixing or elsewhere on the gauge orbit. The 
trick is simply to sample at most once from each 
orbit. Since there are an infinite number of gauge 
orbits (even on the lattice), no finite ensemble 
will ever sample the same orbit twice. This makes 
Gribov copies and gauge-fixing irrelevant in the 
calculation of color-singlet quantities on the lat- 
tice. 

The calculation of gauge- dependent Green's 
functions on the lattice does require that the 
gauge be fixed. The standard choice is naive lat- 
tice Landau gauge, which selects essentially at 
random between the Landau gauge Gribov copies 
for the gauge orbits represented in the ensemble. 
This means that, while the gauge fixing is well- 
defined in that there are no Gribov copies, the 
Landau gauge-fixed configurations are not from 
a single connected FMR. For this reason lattice 
studies of gluon and quark propagators are now 
being extended to Laplacian gauge for compar- 
ison. Laplacian gauge is interesting because it 



is Gribov-copy-free (except on a set of configura- 
tions of measure zero) and it reduces to Landau 
gauge at large momenta. Lattice calculations of 
the Laplacian gauge and Landau gauge quark and 
gluon propagators converge at large momenta and 
hence are consistent with this expectation. 

In conclusion, it should be noted that through- 
out this discussion there has been the implicit as- 
sumption that nonpcrturbative QCD should be 
defined in such a way that each gauge orbit is 
represented only once in the functional integral, 
i.e., that it be defined to have no Gribov copies. 
This is the definition of nonpertubative QCD im- 
plicitly assumed in lattice QCD studies. We have 
seen that this assumption destroys locality and 
the BRS invariance of the theory. An equally 
valid point of view is that locality and BRS 
symmetry are central to the definition of QCD 
and must not be sacrificed in the nonpertuba- 
tive regime, (see, e.g., Ref. j3p|§])- This viewpoint 
impies that Gribov copies are necessarily present, 
that gauge orbits are multiply represented, and 
that the definition of nonperturbative QCD must 
be considered with some care. Since these non- 
pertubative definitions of QCD appear to be dif- 
ferent, establishing which is the one appropriate 
for the description of the physical world is of con- 
siderable importance. 
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